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ABSTRACT 

We use the Least Action Principle to predict the peculiar velocities of PSC-z galaxies 
inside cz = 2000kms~ 1 . Linear theory is used to account for tidal effects to cz = 
15000kms _1 , and we iterate galaxy positions to account for redshift distortions. As the 
Least Action Principle is valid beyond Linear theory, we can predict reliable peculiar 
velocities even for very nearby galaxies (ie cz < 500kms _1 ). These predicted peculiar 
velocities are then compared with the observed velocities of 12 galaxies with Cepheid 
distances. The combination of the PSC-z galaxy survey (with its large sky coverage 
and uniform selection), with the accurate Cepheid distances, makes this comparison 
relatively free from systematic effects. We find that galaxies are good tracers of the 
mass, even at small (< 10 h _1 Mpc) scales; and under the assumption of no biasing, 
0.25 < (3 < 0.75 (at 90% confidence). We use the reliable predicted peculiar velocities 
to estimate the Hubble constant Ho from the local volume without "stepping up" the 
distance ladder, finding a confidence range of 65— 75 kms _1 Mpc _1 (at 90% confidence). 

Key words: galaxies:distances and redshifts - dark matter - distance scale - large- 
scale structure of universe- mcthods:numerical 



1 INTRODUCTION 

The relationship between galaxies and matter is one of the 
most important questions in cosmology. On large scales, 
there is good evidence that galaxies (and in particular IRAS 
galaxies) are good tracers of the mass. For example: there is 
the good alignment of the IRAS and CMB dipoles (Rowan- 
Robinson et al 1990, Strauss et al 1992), the agreement be- 
tween predicted peculiar velocities and Tully-Fisher veloci- 
ties (Willick et al 1997a), and the consistency of POTENT 
density fields with density fields from galaxy surveys (Sigad 
et al 1997). However, on smaller scales, the reliability of 
galaxies as tracers of the mass is far from certain. 

Investigations of peculiar velocities in the Universe have 
relied heavily on use of the linear term in gravitational in- 
stability theory (eg, Rowan-Robinson et al 1990, Strauss et 
al 1992, Willick et al 1997a, da Costa et al 1998). This has 



the advantage of being easy and computationally quick, but 
has the side effect of restricting such analysis to large scales 
~ 10 h~ Mpc, where linear theory is valid. In this paper we 
use the Least Action Principle, which is valid beyond linear 
theory. It is therefore possible to predict peculiar velocities 
for nearby galaxies. The Least Action Principle is similar 
to the more usual action principle, but with the conditions 
that the initial conjugate momenta are zero and the final 
positions fixed, instead of the more usual conditions of both 
the initial and final positions being fixed. It was first used 
by Peebles (1989(P89), 1990(P90)), to study the dynamics 
of the Local Group. Later papers Peebles (1994,1995), pre- 
sented different versions of the method, but still dealt mainly 
with the problem of the Local Group, with only a few ex- 
ternal masses. Shaya et al (1995(S95)), used the method to 
predict motions for galaxies inside cz = 3000kms _1 , and our 
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methodology is very similar to theirs. Recently, Schmoldt & 
Saha (1998), presented a version of the method directly ap- 
plicable to redshift space. 

Previous work on the local velocity field at small scales 
has centred mainly on trying to mimic the observed "cold 
flow" (low value for the rms of galaxy velocities) of local 
galaxies in N-body simulations (Schlegel et al 1994, Gov- 
ernato et al 1997). In both papers it was found to be vir- 
tually impossible for an Qo — 1 CDM model to produce 
a flow as cold of that of the local galaxies. In this paper 
we will compare the observed peculiar velocities of galaxies 
with cepheid distances, to their predicted peculiar velocities 
from the Least Action Principle. This more direct compari- 
son will have the advantage of being less model dependent. 
The emphasis of direct comparisons of predicted and ob- 
served velocity fields has traditionally been on more distant 
galaxies. This is due to the wish to use linear theory and the 
need to have a sufficient number of galaxies to combat the 
large random errors associated with Tully-Fisher distance 
measurements. Our analysis is similar to the one by S95, al- 
though we restrict our attention to a local sample of galax- 
ies (7/12 of our distances are less than lOMpc). We also 
offer two major improvements. Firstly the PSC-z redshift 
survey has the advantages of greater uniformity and larger 
sky coverage, than the optical catalogue of S95. Secondly 
the greater availability of Cepheid distances, gives us a far 
superior (though smaller) set of distances. Both of these im- 
provements should leave our estimate of (3 and other results 
freer of systematic errors. 

It is difficult to measure Ho in the local volume, due to 
the effects of peculiar velocities on the redshifts of galaxies. 
The traditional approach to this problem has been to extend 
the distance scale out to distances where the peculiar veloci- 
ties will be negligible compared to the Hubble flow. However, 
the reliable peculiar velocities which we predict, allow us the 
possibility of measuring Ho locally, without "stepping up" 
the distance ladder. Such a measurement will obviously be 
a useful check on more traditional distance ladder determi- 
nations. 

The layout of this paper is as follows. In the next Sec- 
tion we discuss the Least Action Principle. In Section ^ we 
discuss our mass tracers (mainly the PSC-z galaxy redshift 
survey) and the distance estimations we have used. In Sec- 
tion y we discuss the results from the peculiar velocity com- 
parison. Finally, in Section |^ we give a discussion of our 
results. 



2 THE LEAST ACTION PRINCIPLE 



Xi + 5xi(t), and then integrating equation [j] by parts, 



7 rma -e± + mi v." 



(2) 



To get the equations of motion from the stationary points 
of the action, we require the second term to vanish. Nor- 
mally this is done by setting = at the initial and final 
positions. However, it is also possible to apply the action 
principle with the boundary conditions = 0, at the final 
position and a 2 Sxi/dt — 0, at the initial position. This con- 
dition will be satisfied by the growing mode in linear theory 
and is a natural condition for the early universe. 

We now have to find stationary points of the action. In 
practice, we will only look for minima. The global minimum 
will correspond to the solution closest to the Hubble flow, 
whist other extrema will apply to solutions where (for ex- 
ample) masses have crossed once and our now falling back 
into each other. When Peebles (P90) applied the method to 
the Local Group, only one galaxy (N6822) required a non- 
minimum solution to be well modelled. To find the minima 
numerically we expand the galaxy orbits in terms of a set of 
trial functions, 



Xi(t) = Xf(t ) + ^ C" n fn(t). 



(3) 



Then to find the orbits we minimise the action with re- 
spect to the coefficients Cf >n . The choice of expansion func- 
tions f n {t), is in theory arbitrary; as long as they sat- 
isfy the boundary conditions. Though their growth rate 
should be consistent with linear theory at early epochs (ie 
adx/dt —* a 1112 at a — ► 0). We follow previous papers (eg 
P89, P90, S95), in using the set, 



/»(*) = (! -a)' 



5! 
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1 < n < 5. 



(4) 



The equal spacing of the peaks of these functions allow for 
good convergence and their relatively simple form makes 
them easy to use. The coefficients are found by "walking 
down"; that is by following the gradient of the action to- 



wards the minimum. Normally we start with all C° 
and iterate according to the equation, 

5s 
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Where gf is given by, 

4 G 
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The action S for a system of point masses in an expanding 
universe is (P89,P90,S95), 



/•to fto 

S= dtL = dt 
Jo Jo 



Em»a 2 x? 
2 



- m,i<j)(xi) 



(1) 



Where a is the expansion factor for the Universe (do = 1), 
rm is the mass of the ith particle, x» its position and $>(xj) 
is the gravitational potential. 

The action principle is derived by considering changes 
to the action under infinitesimal changes to the orbits Xj = 



Where p is the mean density of the universe. This process 
con tinues for a sufficient number of iterations (see Section 
4.1), that the Cf„ have converged to near constant values. 
We will then know the orbit of every mass and thus its 
peculiar motion. We place all of the available mass in the 
point masses and unlike previous papers (P90, S95), do not 
consider the possibility of a smooth component to the mass 
(which was found to have little effect on the results). 

So far our discussion has been in terms of point masses. 
In reality, though, we will have to apply the Least Action 
Principle to a set of galaxies. As a result we will have to 
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make the assumption that all of the mass is concentrated in 
these galaxies. An important consequence of this assumption 
is that we r equire the size of galaxy halo (or group halo, see 
Section 3.4) to be less than the average distance between 
galaxies. There is some support for this approximation from 
studies of the Local Group. Zaritsky et al (1989), apply the 
timing argument to the Milky Way- Leo I system and the 
Milky Way- M31 system, and find roughly the same mass 
for the Local Group. As Leo I is at 200kpc and M31 is at 
700fcpc, this implies that the Milky Way does not have a 
large extended halo. In addition, Peebles (1994,1995) finds 
that it is possible to construct satisfactory orbits for Lo- 
cal group satellites and nearby groups with simple compact 
haloes. Another consequence of the assumption that all of 
the mass is concentrated in galaxies is to place limits on the 
distribution of loose mass in the Universe (by loose mass we 
mean mass which is not in haloes). This is particularly im- 
portant in CDM models where up to 50 % of the mass can be 
loose in a fio — 1 universe and greater proportions for lower 
density universes (Governato, private communication). 

At the scale of the Local Group, Branchini & Calberg 
(1994) and Dunn & Laflamme (1995) have tested the Least 
Action Principle in CDM universes. The results from these 
tests are not totally reliable, as the simulations used had 
either been rescaled for use on Local Group, or lacked res- 
olution. Both papers found a tendency for the Least Ac- 
tion Principle to under-predict Qo- Branchini & Calberg at- 
tributed these problems to extended haloes, whereas Dunn 
& Laflamme thought the problem was more to do with loose 
mass. Our analysis is at larger scales, so it is not clear what 
to expect in CDM models. 

The condition that the mass is concentrated at the cen- 
tre of the mass tracers has no chance of being satisfied at 
early epochs. The best we can hope for is that the perturba- 
tions are roughly spherical and centred on the mass tracers. 
We are, though, entirely concerned with the velocities of 
galaxies and these will be generated mainly at late epochs. 
Thus the problems with the tracing of the density field at 
early epochs should not be too serious. 



3 MASS TRACERS AND DISTANCES 

In this Section, we discuss the galaxy catalogues to which 
we apply the Least Action Principle. However, raw galaxy 
catalogues in the form of a list of galaxy positions, redshifts 
and flux, are not suitable to have the Least Action Principle 
directly applied to them. For example, the raw redshifts will 
have a contribution both from the Hubble flow and the pe- 
culiar velocity of a galaxy, so will not be the true measure of 
distance we require. Consequently, this Section will include 
not only a description of the galaxy catalogues we use, but 
also a description of the processes we apply to the galaxy 
catalogues, to make them suitable for calculations with the 
Least Action Principle. 

For peculiar velocity work, there will be certain desir- 
able properties we will want in a galaxy catalogue. The most 
important of these is that the catalogue is uniformly selected 
across the sky, as preferential selection in any region, will re- 
sult in false motions being calculated towards that region. 
Another important property is that we want the catalogue 
to cover as great a fraction of the sky as possible. All cata- 



logues will have an "excluded region" or "mask" around the 
galactic plane. In this region it is impossible to select galax- 
ies reliably. In optically selected surveys this is because of 
absorption from dust. In infrared selected surveys due to 
cirrus emission. It will clearly be desirable for this excluded 
region to be as small as possible. The last desirable property 
is simple density of sampling. The more galaxies we have in 
a given volume, the greater chance we will have of calculat- 
ing accurate peculiar velocities; assuming our galaxies are 
reliable tracers of the mass. 



In Subsection 3.1, we discuss how galaxy bias effects 
our analysis. We then describ e the PSC-z survey, our main 
source of galaxies. Subsection 3.2 is a general introduction 
to the catalogue and is followed by Subsections on: the 
masked regions, the grouping of galaxies, the Virgo Clus- 



ter, tidal effects and redshift distortions. In Subsection 3.8 



we discuss how we have treated local galaxies (those with 
cz < 500kms _1 ). Finally, in Subsection 3.9, we present the 
distance estimations we have used to derive the observed 
peculiar velocities. 

3.1 Galaxy Biasing in our Analysis 

A preliminary issue to be discussed is the effect of bias in 
the galaxy distribution on our analysis. We need a catalogue 
of mass tracers on which to apply the Least Action Princi- 
ple, this will be a catalogue of galaxies and galaxy groups. 
This means that our analysis will depend on the relation- 
ship between the galaxies and the mass. In addition to our 
requirement that the mass is concentrated in the galaxies 
(see Section^), we have uncertainty from whether the mass 
associated with galaxies is dependant on their environment. 
This is the problem of the bias of the galaxy distribution. 
The standard approach to this problem is to use the linear 
bias factor b, to quantify the uncertainty in this relationship 
(b is given by 5 g = b8 m , where S m and S g are the over- or 
under-densities in the galaxies and the mass respectively). 
This has the advantage that in linear theory, the effect of b 
is degenerate with that of Qo and the results only depend on 
(3 — Q l Q 6 /b. However, the Least Action Principle does not 
rely on linear theory and one of the main reasons for using 
it is to probe small scales at which linear theory is invalid. 
Thus we can no longer expect Qo and b to be fully degen- 
erate. We chose to weight our galaxies on the assumption 
that there is no bias (ie 6=1). This is not as restrictive as 
it might appear. Although strictly Qo and b are not degen- 
erate, we expect that even at non-linear scales a large de- 
generacy will persist. Support from this comes from (Dekel 
et al 1993), who investigated the peculiar velocity field us- 
ing the POTENT methodology. This methodology included 
non-linear effects and could in principle break the degener- 
acy between Qo and b. They found, though, the individual 
constraints found on fio and 6 were extremely weak and only 
the parameter j3 was constrained to good accuracy. 

The POTENT analysis is at fairly large scales, with the 
galaxy and velocity fields being smoothed on 1200 kms -1 . 
Our analysis is at smaller scales and is thus more likely to 
see the degeneracy between Qo and b broken. Consequently, 
we have investigated the effect of bias in our calculations. 
We can do this by modifying the expression used for the 
potential in equation |l]. Including bias, the potential can be 
expressed as 
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Pmb 



-Ga 2 



Pb 



(8) 



Where p(x') is the density of galaxies at x' and pb is the back- 
ground density of galaxies. p m b is the background density of 
matter and the factor (p m b/pg) just introduces a physical 
mass (as oppossed to galaxy number) into the equation. For 
b=l, this reduces to the standard form for the potential 
(Peebles 1980). Using equation | for the potential, equation 
W becomes, 



-TvGpax.', 



■ x, 



(9) 



We can then use equation ^ in equation ^ and find the min- 
ima of the action for a biased galaxy distribution and thus 
investigate the effect of bias. We should note, though, that 
our treatment is not entirely satisfactory. We have intro- 
duced bias in equation | in a manner that is reliant on 
smooth and continuous density fields. However, the Least 
Action Principle applies to a set of point masses not contin- 
uous fields. We have thus introduced bias in an inconsistent 
manner. To have introduced bias consistently would of re- 
quired us to weight galaxies according local density. Such a 
process would be both complicated and noisy and thus we 
have chosen not to pursue it. We can, though, expect an 
investigation of bias using equations ^ and ^ to give at least 
an indication of its effect on our results. 

We thus did a series of runs with one of our input cata- 
logues with different fio and b but constant j3 (the catalogue 
we used had a PS C-z local weighting and a cloned mask see 
Subsections 3.2 - 3.8). We then compared the magnitudes 
of the predicted velocities in the different calculations. For 
/3 = 0.5 (fio = 0.3 b = 1) then reducing b to 0.5 (fi = 0.1) 
resulted in the velocities decreasing on average to a factor 
0.77 of the b — 1 velocities. Whilst, increasing b to 2 in- 
creased the velocities to a factor 1.29 of the 6=1 velocities. 
Given these results, it is clear that we still have a large de- 
generacy between fio and b and the correct quantity for us 
to estimate is (3 . Indeed, with the fairl y we ak limits we will 
eventually obtain on j3 (see Subsection 4.3) it would be fu- 
tile to attempt to estimate fio and b independently. Thus 
we will only estimate /3 and not fio and b. We will, though, 
quote are initial results in terms of fio as we do not wish 
to hide the assumption that 6 = 1, used in obtaining these 
results. 

An additional consideration comes from the possibilty 
of non-linear bias. If non-linear terms are present in the bias- 
ing relation, then they are predicted to be more important 
at smaller scales (Mo & White 1996) and so could poten- 
tially be important for our work. However, most non-linear 
biasing models, tend to differ substantially from linear bias 
only in clusters or voids (see the models of Mann, Peacock & 
Heavens 1998) . The local volume is neither of these (Kraan- 
Korteweg & Tammann 1979), so for our analysis most non- 
linear biasing schemes will be consistent with linear bias and 
our weighting scheme. As a result our conclusions on how 
galaxies trace the mass, will be for a general trend and not 
a particular biasing scheme. 



3.2 The PSC-z Survey 

The galaxy catalogue we use is the PSC-z survey (Saun- 
ders et al 1996). This is a redshift survey to the full depth 
(Seo > 0.6 Jy) of the IRAS Point Source Catalogue. It has all 
the properties required for accurate calculation of peculiar 
velocities. Selection from the IRAS Point Source Catalogue 
ensures near uniform selection across the sky. In particu- 
lar, it avoids the problems with optically selected surveys of 
matching the northern and southern skies. The PSC-z sur- 
vey also has excellent sky coverage, with 84.1% of the sky 
in the survey. The flux limit of 0.6Jy makes the PSC-z the 
deepest complete all sky IRAS redshift survey, ensuring that 
we have sufficient density of sampling to accurately calculate 
peculiar velocities. 

The Least Action Principle is computationally ex- 
pensive to apply, so we only use galaxies with cz < 
2000kms _1 (unless otherwise stated all of our redshifts have 
been transformed to the Local Group frame using the Yahil, 
Tammann & Sandage et al (1977) correction). For an Ho 
of 80 kms _1 Mpc _1 , our furthest distance of 17.8Mpc cor- 
responds to cz — 1420kms _1 , so we have a reasonable 
buffer zone to the edge effects at cz = 2000kms -1 . Galax- 
ies with cz < 500kms _1 , or within 12° of M87 (ie in the 
Vir go cluster) are treated separately (see Sections 3.S and 



3.5 respectively). We impose a luminosity cut, to avoid any 
weighting problems with the poorly determined faint end 
of the luminosity function. The luminosity cut is such that 
all galaxies would have a flux of 0.6Jy or greater when 
placed at cz = 1060kms~ 1 . We give each galaxy a weight 
w{cz) = l/tp(cz), where 4>(cz) is the PSC-z selection func- 
tion (Saunders et al 1999). The selection function quantifies 
the loss of galaxies with distance, due to the flux cut. Its 
value is the number density of galaxies that would be ob- 
served at a given redshift in the absence of galaxy fluctu- 
ations. The number of galaxies in this main sample is 929. 
This includes galaxies below the luminosity limit, as these 
may become sufficiently luminous when we iterate t he p osi- 
tions to account for redshift distortions (see Section 3.7). 



3.3 Treatment of Masked Region 

Although the PSC-z mask is small, we still have to account 
for it. We do this by filling the sky behind the mask with 
fictitious galaxies. To get a feel for the errors involved, we 
adopt two very different schemes for doing this. The first is 
the "random mask" , where galaxies are placed at random in 
the masked region, with a space density consistent with the 
PSC-z selection function. The second is the "cloned mask" 
where we extrapolate structures across the Galactic plane, 
from the galaxies we observe above and below the plane. 
Masked regions away from the Galactic plane (mainly the 
small areas where IRAS did not survey) are still filled with 
random galaxies. For the random mask we add a total of 169 
"galaxies" and for the cloned mask we add 145 "galaxies" . 



3.4 Grouping of the Galaxies 

As we described in Section we are looking for minimum 
solutions to equation ^[ as the global minimum is the solu- 
tion closest to the Hubble flow. This solution is likely to be 
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adequate, except in regions at the centres of groups and clus- 
ters. In these regions, the motions of galaxies are likely to 
be a general stationary point of the action, as the orbits will 
have crossed (see P89, P90). So to avoid problems with these 
stationary solutions, we group our galaxies and just follow 
the motion of the centre of mass. The grouping of galaxies 
also helps us to deal with redshift distortions. We attempt 
to iterate out the coherent component to the redshift dis- 
tortions during our calculations (see Section 3.7). However 



tions using linear theory and add them to the Least Action 
Principle results. We use the standard equation, 



there is nothing we can do about the incoherent component 
to redshift distortions generated on small scales. By group- 
ing the galaxies we can hope to average out the worst effects 
of these incoherent redshift distortions. To group we use a 
"Friends of Friends" algorithm, similar to Huchra & Geller 
(1982). The linking is Ad = 0.6 h _1 Mpc tangentially, and 
Az — 150kms _1 radially, for galaxies with average redshift 
inside cz — 1060kms~ 1 . For galaxies outside the volume lim- 
ited sample we scale the lengths, to allow for the decreasing 
density of selected galaxies. This is done according to, 



Ad = 0.6 ( *g™> Yh-^Mpc, 
Az = iBof^Vfcm.- 1 . 

V n cz > J 



(10) 

(ii) 



The radial linking is so much greater to allow for the effects 
of velocity dispersion inside the groups. Once a galaxy is 
assigned to a group, it is considered to be at the centre 
of mass of that group. We have compared the groups with 
those of the Tully (1987) catalogue. Most of the time the 
two agree fairly well, though there were differences, but we 
are confident that the grouping is working well. 



3.5 The Virgo Cluster 

The Virgo Cluster is the most important object for the dy- 
namics of the local volume, so it is important to treat it well 
in our calculations. To do this, we have treated it separately 
from the rest of the galaxies. This is necessitated by the 
large size and velocity dispersion of the cluster. The link- 
ing lengths of the grouping algorithm, have been chosen for 
small to medium groups, and are thus inappropriate for the 
Virgo cluster. We also need to take account of blueshifted 
galaxies, which outside of the Local Group, are only present 
in the Virgo Cluster. We consider all galaxies within 12° of 
the Virgo core (which we take to be M87), to be a part of 
the Virgo cluster. Provided that they have a redshift below 
cz = 2000kms _1 , galaxies with blueshifts in this region of 
the sky are also placed in the Virgo cluster. We cannot ex- 
tend the upper redshift limit beyond cz — 2000kms _1 due to 
contamination from background galaxies. On these criteria, 
105 galaxies are placed in the Virgo Cluster, which we place 
at the angular position of M87. For its redshift we use the 
value given by Bingelli et al (1987), not including the dwarf 
galaxies. 



3.6 Tidal Effects 

Galaxies outside cz = 2000kms~ 1 , could potentially gener- 
ate large tidal motions. We thus calculate these tidal mo- 



V = — — EiWi^-. 
47r r J 



(12) 



Where iu, is the weight of the ith galaxy (see Section 3.2). 
The tidal velocities are calculated with all galaxies between 
cz = 2000kms" 1 to cz = 15000kms _1 in the PSC-z. The 
galaxies are placed in the Local Group frame if they are 
inside cz — 3000kms _1 , but in the CMB frame if they are 
outside. We use a random mask model throughout for the 
tidal effects. The tidal velocities are also added to the pecu- 
liar velocities when we compensate for redshift distortions. 
Clearly the positions of our galaxies, with just the switch 
in frames at cz — 3000kms _1 are very crude. We have com- 
pared the tidal velocities calculated using these positions, 
with those obtained using a full reconstruction of real space 
positions. The results compare very well, with most veloci- 
ties differing by less than 10kms~ . 



3.7 Redshift Distortions 

The redshifts of the galaxies and groups will suffer from red- 
shift distortions due to their peculiar motions. The simple 
use of these redshifts as distances would cause large errors. 
As a result we attempt to iterate out the effects of redshift 
distortions in a manner similar to S95. This involves iter- 
ating the position of a galaxy/group so that the quantity 
l/c(H d + v), is equal to the redshift of the galaxy (where 
d is the distance to the galaxy and v its radial peculiar ve- 
locity, relative to the Local Group). The peculiar velocities 
are calculated using the Least Action Principle (see Section 



1). We iterate the positions 10 times, and the distance a 
galaxy is allowed to move is slowly scaled down throughout 
the 10 iterations. By the last iteration most of the residues 
are small ( < 0.5 h -1 Mpc in most cases). 

We have proper distances for most of the local galaxies 
(ie those with cz < 500kms _1 ). So, in order to avoid the va- 
garies of iterating out the redshift distortions we place these 
galaxies, where we can, at our distance estimates. It is very 
important to place these galaxies at the correct distances, 
as due to their close proximity, they can potentially have a 
large effect on the motion of the Local Group. The motion 
of the Local Group is crucial as all velocities are measured 
relative to it. To be consistent we also place the galaxies for 
which we have Cepheid distances outside cz = 500 kms - , 
at those Cepheid distances. An exception is Virgo, which 
is allowed to iterate to near its distance, and then placed 
at its distance. This is to allow the correct iteration of the 
redshifts of galaxies near the Virgo Cluster. Although the 
grouping is done with the galaxies in the Local Group frame, 
we then attach galaxies below the luminosity cut to groups, 
or treat them as massless particles. Then, as we iterate out 
the redshift distortions, and change the distance to galax- 
ies, we re-assess whether galaxies are above or below the 
luminosity limit and act accordingly. This may seem quite a 
trivial effect, but a substantial number of galaxies do change 
status, in a very systematic way across the sky. 
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Table 1. PSC-z and Optical weightings for the local galaxies. 
The total weight in each column is 1.00. The group names refer 
to those of Kraan-Korteweg & Tammann (1979). We give alter- 
natives names in brackets, normally the brightest galaxy. 



3.8 Weighting of Local Galaxies 

The galaxies inside cz = 500kms _1 are treated differently 
from the rest of the survey. This is partly to reflect our in- 
creased knowledge of these galaxies and partly due to their 
importance in determining the Local Group's velocity (rel- 
ative to which all other velocities are observed). We did not 
apply our grouping algorithm to these galaxies, but instead 
use the 7 groups denned by Kraan-Korteweg & Tammann 
(1979), but follow the pleading of Schmidt & Boiler (1992b) 
for a subdivision of the B7 group. 

We have two distinct ways of assigning weights to these 
local galaxies. The first is using the PSC-z catalogue, each 
group or galaxy is given a weight according to how many 
galaxies it contains above the luminosity limit; when placed 
at its distance for Ho = 80kms _1 Mpc _1 . The actual value 
of Ho used makes little difference to the weights. If a galaxy 
has no distance and is not in a group with a distance, we use 
its redshift. A second weighting is given by optical luminos- 
ity. We take the optical luminosity data from the catalogue 
of local galaxies compiled by Schmidt & Boiler (1992a), each 
galaxy and group is weighted in proportion to its luminos- 
ity and and all galaxies and groups with a luminosity above 
lO lo Z/0 are included, (B7a is just under this limit, but too 
close to the Local Group to cut). To tie these relative weights 
into the PSC-z weightings used for the rest of the galaxies 
we multiply the relative weights by the total weight of PSC- 
z galaxies inside cz = 500kms _1 for Ho = 80 kms _1 Mpc _1 . 
The results of these schemes can be seen in Table [jj. The 
main difference in the two weighting schemes, is that while 
the PSC-z weighting tends to be fairly even across the galax- 
ies and groups, the optical weighting tends to favour a few 
large groups. With the two possible mask compensations 
and the two possible local weightings, we have four possible 
catalogues. 

We have two main aims in having two weighting 
schemes. Firstly, we wanted to choose the weighting scheme 



which optimised the results, and secondly, to see how robust 
the results are to different weighting schemes. The issue of 
robustness is important. Although either of these weight- 
ing schemes might be on average a good tracer of the total 
mass, we expect that the weight assigned to any one galaxy 
or group is subject to large errors. We can gauge the pos- 
sible effects of these errors by looking at the results of the 
two weighting schemes. 



3.9 Distances 

We present our distances in Table []. All the distances are 
derived from Cepheids and seven of the 12 are inside lOMpc, 
making our analysis sensitive to small scales. All of the 
galaxies are in the PSC-z catalogue. The errors consist of 
the quoted distance errors given in the references added in 
quadrature to a "projection" error for groups. The projec- 
tion error accounts for the fact that the galaxy to which we 
have a distance may not be at the centre of the mass of the 
group. This error is given by 



Ad: 



V2n 



(13) 



Where the AOi are the angular distances of the galaxies 
from the centre of mass, n is the number of galaxies in the 
group and d is the distance to the group. This error is nor- 
mally the minor partner, only dominating for B5, B6 and 
B7a. It is possible that our distances are affected by system- 
atic errors. The most likely sources are either problems with 
the calibration due to errors in the distance to the LMC, or 
errors due to metallicity effects. Problems with the calibra- 
tion at the LMC will affect all the distances by the same 
fractional amount. Thus, although it will have direct effect 
on the estimate of Ho, it will leave the rest of our conclu- 
sions unchanged. Significant metallicity effects would have 
the potential to effect all of our results, but luckily they are 
thought to be small (Kennicutt et al 1998). 

To allow us to better combat redshift distortions in the 
local volume, we place the groups Bl, B4 and B7b and the 
galaxy N3115 at non-Cepheid distance measures. To avoid 
possible systematic errors these are not used in the pecu- 
liar velocity comparison. Bl is placed at 4.2Mpc, from the 
K-Band surface brightness fluctuations measured by Lup- 
pino & Tonry (1993). Graham et al (1982) places B7b IMpc 
more distant than B7a, so we place it at S.lMpc. The Tully- 
Fisher ratios of B4 and B5 from Willick et al (1997b), led 
us to place B4 at 3.7Mpc. N3115 is placed at 10.8Afpc from 
the Planetary Nebula Luminosity Function determination of 
Ciardullo et al (1993). 



4 RESULTS 

In this Section we discuss the results from our velocity com- 



parison. The first Section 4.1, deals with the exact details of 
how we apply the Least Acti on P rinciple. We then discuss 
the goodness of fit in Section 4.2 and the confidence limits 
we obtain on Qq and Ho in Section 4.S . 
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Table 2. A Table of distances to nearby galaxies and groups, the 
redshifts are in the Local Group frame with the Yahil et al (1977) 
correction. 
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ms 1 
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6% 
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100% 
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100% 
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97% 


98% 
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Table 3. The convergence of velocities with number of iterations. 
The velocities are compared to those obtained with 1000 itera- 
tions. The results here are for Co = 0-3, Hq = 70kms~ 1 Mpc~ 1 , 
PSC-z local weighting and cloned mask. Results from other mod- 
els are very similar. 



4.1 Application of Least Action Principle 

As explained in Section ^, we iterate the Cf >n towards min- 
imum using equation ^. A crucial point is how many iter- 
ations are required to achieve good convergence. We have 
approached this problem in a different manner to previous 
papers (P90, S95), who looked at the changing value of the 
gradient and stopped calculations when this was sufficiently 
small. Instead we compared the velocities predicted after a 
number of iterations, with those obtained after 1000 itera- 
tions. Typical results can be seen in Table [j. As the table 
shows, initial convergence to rough values is very quick, but 
final convergence takes a great many iterations; with veloc- 
ities only changing by a small amount each iteration. We 
thus choose to iterate the C™ n 250 times, iterating galaxy 
positions every 10 C" n iterations for the first 100 iterations 
(see Section 3.7). The rough velocities used to iterate the 



positions are easily adequate; as we only aim to drive the 
residues to < 50fcms _1 . 

In integrating equation ^, we use 10 timesteps. Compar- 
isons with results calculated with 40 and 100 timesteps, sug- 
gests that this gives us ~ 85% of velocities within 20kms~ 
and ~ 98% of velocities within 50fcms _1 . As convergence is 
slightly slower for a larger number of timesteps, it is clearly 
preferable to have (as we do) a low number of timesteps com- 
bined with a large number of iterations, (the computational 
time is directly proportional to the number of timesteps mul- 
tiplied by the number of iterations). 

As mentioned in Section ^, the solution for which we 



seek is the global minimum. This should give us the galaxy 
orbits which are closest to the Hubble flow, and is the nat- 
ural extension of linear theory. The existence of a global 
minimum is guaranteed as S — > oo as C"„ — > oo. However, 
it is also possible that other local minimum exist. Indeed 
trials from large random values of C" n reveal that multiple 
minima exist. As a result approximately 15 % of galaxies and 
groups can have two distinct solutions (by distinct solution 
we mean a difference in velocities of greater than 50 kms -1 ). 
Unfortunately, it is impossible to calculate the differences in 
the value of the action between solutions accurately enough 
to ascertain which solution corresponds to the global mini- 
mum. For pure Hubble flow we would have C" n — 0, so it 
seems reasonable that the solution obtained when we iterate 
from C" n — is the global minimum. Certainly, the use of 
this solution has become conventional in Least Action work 
(ie S95). T hus it is this solution which we use for the results 
in Sections h. 2] and 4.3. A particular problem we experience 



is that the Local Group has two distinct solutions for three 
of the catalogues. Only for the optical local weighting and 
random mask catalogue do the trials from large random val- 
ues of C" n fail to reveal a second solution. For the other 
three catalogues, the two Local Group solutions have the 
potential to cause great problems as all velocities are mea- 
sured relative to the Local Group. We have, by fixing the 
initial values for C" n , explored the other Local Group solu- 
tions for the three catalogues. The results found with these 
solutions are consistent with those presented below, with 
similar confidence limits being obtained. For example with 
the PSC-z local weighting and cloned mask we obtain limits 
0.1 < fi < 0.40 and 65 < H < 72.5 kms^Mpc" 1 with the 
"standard" solution, iterating from C" n = and limits of 
0.05 < Q < 0.25 and 65 < H < 70 kms^Mpc -1 with the 
"second" solution. 



4.2 Goodness of Fit 

A qualitative impression of the goodness of fit can be ob- 
tained from Figure Q This presents scatterplot diagrams of 
the predicted velocities from the Least Action Principle cal- 
culations against the observed velocities from the Cepheid 
distances. All of the plots show results from calculations with 
a PSC-z local weighting and Cloned mask. The error bars 
are given by a i = HqAD. Where AD is the error in the dis- 
tance and a i is the resulting error in the observed velocity. 
For all our results we have calculated the observed veloc- 
ity errors using Ho = 70fcms _1 Mpc _1 ; we take no account 
of errors in the predicted velocities. For all the scatterplots 
where Qo and Ho are near best fit, (ie around Qo = 0.2, 
Ho — 70kms~ 1 Mpc -1 , see below) almost all of the points 
lie on the observed velocity equals predicted velocity line, to 
within their errors. The only exception is the 7V1023 group 
(marked in the bottom left plot), even at best fit this point is 
a considerable distance from the line. So with the exception 
of this galaxy we are getting a very good fit. 

To make our results more quantitative we performed a 
reduced x 2 test, 

2 _ 1 — u< ") 2 



N- M ^ 



(14) 



Where v P i is the predicted velocity from the Least Action 
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Figure 1. Scattcrplots showing the effect of varying/fo and Qo- All results are with PSC-z local weighting and a cloned mask model. 
The line shown is observed velocity= predicted velocity. 



Principle, v i is the observed velocity and a a i its error. N 
is the number of data points and M the number of fitted 
parameters. 

Predicted velocities using the Least Action Principle 
where obtained on a grid of f2rj and Ho for the four possible 
input catalogues. In Table ^ we present fio, Ho and the value 
of the reduced x 2 statistic at best fit for the four catalogues. 
We do this using all the galaxies and excluding the outlier 
N1023 group. 

For the reduced x 2 test we are performing the 95% con- 
fidence limit is 1.83. Thus the results confirm the impression 
of the scatterplots that the fit is very good, especially when 
the AT023 group is discarded. The PSC-z local weight is 
performing better than the optical weight and by a smaller 



margin the cloned mask is preferred to the random mask. 
Importantly we see that the results are reasonably robust to 
both local weighting and mask compensation, with the best 
fit Oo and Ho being relatively close for the four models. 

We are not sure why we are having such problems with 
the AT023 group. It is not in the 15% of galaxies and 
groups for which trials from large random values of C" n 
reveal multiple solutions. Nor is the predicted peculiar ve- 
locity strongly effected by the multiple solutions of the Local 
Group. It lies reasonable close to the galactic plane, so we 
could be experiencing problems in how we compensate for 
the mask. The group is an outlier for both the mask mod- 
els, but it is possible that both models miss an important 
structure (ie a reasonable compact group) that is essential to 
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Table 4. Best Fit x 2 values, with and without N1023. The values of Oq and Hq given are those at the best fit \ 2 ■ 



accurately predict the velocity of A1023. Another possibil- 
ity is a problem with the distance to the group. This seems 
less likely as the Cepheid distance agrees well with other 
measurements of 9.6Mpc from the Surface Brightness Fluc- 
tuation method (Tonry et al 1997), and W.lMpc from the 
Planetary Nebula Luminosity Function method (Ciardullo 
et al 1991). For the purposes of the next Section, where 
we calculate confidence intervals, we chose to exclude the 
iV1023 group. If either of the problems mentioned above are 
true, then clearly we are well justified in doing so. However, 
even if JV1023 is a natural outlier then it is necessary to ex- 
clude it. The small number of galaxies that we have are not 
robust to its effect and fairer limits are obtained without it. 



4.3 Confidence Intervals 

The velocities that we predict for galaxies and groups de- 
pend mainly an the input value of fio- There is also a weak 
dependance on Ho as this changes where we place the galax- 
ies with measured distances. The main effect of Ho in the 
peculiar velocity comparison, is that it changes the observed 
velocities (given by v — cz — Hod) not the predicted veloc- 
ities. Figure ^ shows the effects changing of fio and Ho- 
Changing fio stretches or compresses the range of predicted 
velocities, thus for low fio our range of predicted velocities 
is too narrow, whilst for higher Qo it is too broad. Chang- 
ing Ho shifts the observed velocities up and down, thus for 
low Ho the observed velocities tend to be too large, whilst 
for high Ho they tend to be small. In effect we can think of 
fio as changing the gradient of our data and Ho as chang- 
ing the zero-point, only for the certain values of Qo and Ho 
will our data be consistent with the observed velocity equals 
predicted velocity line. 

We can make a quantitative assessment of the 90 % 
confidence limits by using a Ax 2 test. The results of this is 
shown in Table ^. For all four models, we get similar ranges 
of Qo and Ho . This indicates that we are robust to the Local 
weighting scheme and to the errors we expect in weighting 
given to any individual galaxy or group. We are also rea- 
sonable robust to mask compensation model, though the 
random mask does allow for a higher upper limit on fio. 
We take our final confidence limits from the union of the 
results with the PSC-z local weighting, as this weighting 
was preferred. This gives 65 < Ho < 75 kms Mpc , and 
0.1 < Qo < 0.60. We can convert to /3 = Qo' 6 /b, under the 
assumption that 6 = 1 (see Subsection 3.1). This gives us 



0.25 < /3 < 0.75. As in the determination of Ho the distance 
to Virgo is somewhat controversial we have calculated an Ho 
confidence interval without using this distance for the PSC-z 
local weighting and cloned mask catalogue. We get the same 
range as before (ie 65 < H < 72.5 kms" 1 Mpc" 1 ), so our 
estimate does not depend very heavily on the distance to 



Mask 
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random 

cloned 

random 

cloned 



Local 

Weightings 

PSC-z 

PSC-z 

optical 

optical 



Ho Range 

klTlS 1 XI jir 

67.5- 75.0 
65.0- 72.5 
65.0- 72.5 
65.0- 72.5 



Co range 

0.10- 0.60 
0.10- 0.40 
0.05- 0.40 
0.05- 0.30 



Table 5. The preferred ranges of Qq an< i H) 



Virgo. For interest, we tend to calculate an infall into Virgo 
of around 200kms" 1 , for p = 0.5. 



5 DISCUSSION 

Using the Least Action Principle and a catalogue of galaxies 
drawn from the IRAS PSC-z survey, we have predicted re- 
liable velocities for nearby galaxies. This is consistent with 
the idea that PSC-z galaxies are good tracers of the mass at 
small as well as large scales. As mentioned in Section |^, we 
are assuming galaxies are point masses and would not have 
expected to predict reliable velocities if galaxy haloes where 
very extended, or the universe contained a lot of loose mat- 
ter. Consequently our ability to predict reliable velocities on 
scales ranging from 2Mpc to 15Mpc, argues against galaxies 
having extended haloes beyond 2 — 3Mpc. It is also difficult 
to reconcile large amounts of loose matter with the good 
fit we observe. A possible solution to this is that the loose 
matter is distributed fairly evenly so that it produces no 
net accelerations. Another possibility is that the clustering 
(above galaxy scale) of the loose matter is tightly correlated 
with the clustering of galaxies. CDM models contain a lot 
of loose matter (see Section |^) , so the good fit we observe 
is a challenge for such models. Indeed as discussed in the 
introduction the two previous attempts at testing the LAP 
in CDM simulations (Branchini & Calberg 1994, Dunn & 
Laflamme 1995) were not successful. It would clearly be de- 
sirable to repeat these investigations with an increased num- 
ber of particles and consequently higher resolution. It could 
then be established whether the LAP should predict reliable 
velocities in CDM universes and allow an assessment of the 
effect on LAP predictions of galaxy merging, mass transfer 
and tidal disruption. 

If we neglect bias (ie assume 6=1) then the Qo results 
confirm the difficulty of fitting the local flow for values of 
Qo around unity. Governato et al (1997) found no cases of 
a CDM simulation with either Qo — 0.3 or fio = 1, produc- 
ing a "cold flow" as low as that observed. Their definition 
of "cold flow" was the rms velocity inside cz = 500kms _1 . 
They used a value of 60kms~ 1 from Schlegel et al (1994). 
The data we use here gives an equivalent "cold flow" of 
43 ± 35kms _1 , so the results are consistent. However our 
results have greater model independence, especially with re- 
gards to the normalisation of the power spectrum. 
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Direct comparisons of peculiar velocities have been 
made at larger scales, with the use of Tully-Fisher distances. 
Willick et al (1997a) use 1.2Jy galaxy redshift survey and 
the Mark III catalogue of peculiar velocities to achieve an 
estimate of /3 = 0.49 ± 0.07. da Costa et al (1998) use the 
1.2 Jy survey with SFI peculiar velocities to get (3 = 0.6±0.1. 
The confidence limits we derive of 0.25 < < 0.75, are 
quite consistent with these estimates. For the PSC-z local 
weighting we get a best fit Qo of around 0.3, which gives 
f3 = 0.5. This agrees well with Willick and da Costa es- 
timates and suggests that there is little variation of the 
measured (3 with the scale the estimate is made on. Such 
a dependence of measured with scale has been invoked 
to explain the difference in (3 estimates between velocity- 
velocity comparisons (Willick et al 1997a, da Costa et al 
1998, and this paper) with density-density comparisons us- 
ing the POTENT methodology (Sigad et al 1998, they es- 
timate f3 = 0.89 ± 0.12). As we find no evidence for such 
a scale dependence, we suggest the source of the discrep- 
ancy must lie elsewhere. S95 use the Least Action Principle 
and Tully-Fisher data to estimate flo = 0.17 ± 0.1 (giving 
(3 — 0.35 ± 0.15, for b = 1). Again this is consistent with 
our results, but despite the similarities in methodology, our 
results are less comparable with this paper, which uses an 
optical catalogue for its mass tracers. 

The classical way to estimate Ho, has been through 
the distance ladder. Typically this means using Cepheid dis- 
tances to nearby galaxies, groups and clusters to calibrate 
either type la supernova (SN la), or the Tully-Fisher rela- 
tion. Direct estimates at the Virgo cluster are also made. 
There tends to be a bifurcation of estimates, for example 
Freedman et al (1998) give a value Ho — 73 ± 6 (statisti- 
cal) ±8 (systematic) kms _1 Mpc~ , whilst Tammann (1998) 
gives Ho = 56 ± 7kms~ 1 Mpc~ 1 . The disagreements caus- 
ing these different estimates are wide-ranging. Some of the 
more important are: the distance to the Virgo cluster, the 
nature of the Fornax cluster, suitable SN la observations for 
calibration purposes, the use of the decline rate- absolute 
magnitude relation for SN la and the magnitude of system- 
atic effects when using the Tully-Fisher relation. We obtain 
the confidence interval 65 — 75 kms _1 Mpc _1 , which favours 
the higher estimate. Importantly the methodology we have 
used by-passes many of the disagreements, relying mainly 
on the noncontroversial Cepheid distances. Of the disagree- 
ments listed above only the distanc e to Virgo applies to our 
estimate, and as shown in Section 4.3, it is of little impor- 
tance to our estimate. 

Recently a number of Ho estimates from the HST 
Key Project have been presented. These Key Project es- 
timates have been made with a variety of secondary dis- 
tance indicators. These include: the Surface Brightness Fluc- 
tuations method (Ho = 69 ± 4 (statistical) ±5 (system- 
atic) kms _1 Mpc~ 1 , Ferraresa et al 1999), the Fundamen- 
tal Plane relation Ho = 78 ± 7 (statistical) ±8 (systematic) 
kms^Mpc" 1 , Kelson et al 1999), the Tully-Fisher relation 
(Ho = 71 ± 4 (statistical) ±10 (systematic) kms _1 Mpc _1 , 
Sakai et al 1999) and the Type IA SN method H = 68 ± 2 
(statistical) ±5 (systematic) kms~ 1 Mpc _1 , Gibson et al 
1999). Perhaps the most interesting of these estimates is 
the Type la SN result (Gibson et al 1999). This method has 
traditionally favoured the lower Ho estimates (e.g. Saha et 
al 1997,1999) and thus this higher estimate removes one of 



the major objections to a high value for Ho- Our estimate 
for Ho (65 — 75 kms _1 Mpc _1 ) agrees very well with these 
HST Key Project determinations. Given that our estimate 
avoids many of the traditional disagreements in distance lad- 
der determinations and the HST Key Project aims to pro- 
vide a definitive estimate of Ho, this amounts to a strong 
argument for a high Ho. 

Although our estimate is made in the local volume, it 
is not a local estimate of Ho- We use the PSC-z selection 
function to weight the galaxies, so if we live in a local over- 
or under-density this will be fully taken into account. As a 
result any net expansion or contraction of the local volume 
will be accounted for in the Least Action Predictions. Thus 
our estimate of Ho will in effect be made to at least the aver- 
gage depth of the PSC-z (cz — 8100 kms -1 ) and will be a 
global not local estimate. We would like to reiterate that like 
all estimates of Ho made with Cepheids, this value is vulner- 
able to errors from metallicity effects, or more importantly, 
an error in the distance to the LMC. 
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